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摘 要 : 基于 中 国 第 二 次 冰川 编目 数据 、 国 家 青藏 高 原 科学 数据 中 心 的 降水 数据 和 气象 站 的 实测 气温 数据 并 结合 
数字 高 程 模型 ,利用 加 入 辐射 项 的 分 布 式 度 日 模型 ,以 已 有 监测 冰川 资料 为 基础 获取 度 日 因子 ,通过 优化 参数 ,对 
祁连山 流域 尺度 1961 一 2013 年 的 冰川 物质 平衡 变化 序列 进行 重建 。 绪 果 表 明 :祁连山 的 冰川 物质 平衡 变化 以 1995 
年 为 界 可 以 划分 为 2 个 阶段 。1961 一 1995 年 ,祁连山 冰川 物质 平衡 正 负 交替 变化 日 变化 幅度 较 小 ,冰川 的 消融 量 和 


积累 量 基本 持平 ， 


Is 


平均 物质 平衡 为 (0.11+0.13)m w.e. +a”; 1996—2013 4E , VK JI TEE Filios op , YO 8 RY 


物质 平衡 为 (-0.54+0.13 )m w.e…a ,主要 原因 是 正 积 温 的 逐年 升 高 导致 冰川 消融 加 剧 。1961 一 2013 年 祁连山 冰川 


的 累积 物质 平衡 为 (-12.76+4.24)m w.e. ,空间 分 布 上 祁连山 


关键 词 : 祁连山 ; 流域 ; 冰川 物质 平衡 ; 重建 


冰川 物质 平衡 作为 气候 的 产物 ,其 动态 变化 会 
引起 冰川 末端 位 置 、 面 积 、 冰 储量 及 径流 的 变化 "。 
随 着 气候 变 暖 ,全球 大 多 数 的 冰川 呈现 退缩 状态 ， 
并 且 在 最 近 20 a 出 现 了 加 速 的 趋势 “' 。 研 究 表明 
在 1.5 % 气 候 情 景 下 ,到 21 世纪 末 亚 洲 高 山区 的 冰 
川 储量 只 有 目前 的 64%+7%"。 山 地 冰川 作为 一 种 
淡水 资源 ,在 农业 用 水 .水力 水 生生 物 以 及 基本 的 
水 供应 方面 起 着 重要 的 作用 *”。 尤 其 是 在 干旱 半 
干旱 区 ,夏季 冰川 融 水 是 河流 的 主要 补给 来 源 ,对 
径流 具有 调节 作用 。 

祁连山 的 冰川 区 是 某 达 木 盆地 干旱 半 干 时 区 
与 河西 走廊 的 “ 湿 岛 ,监测 和 估算 祁连山 区 冰川 
物质 平衡 的 变化 对 干旱 区 的 意义 重大 。 冰 川 物质 
平衡 的 重建 能 够 恢复 其 长 时 间 序 列 的 变化 ,有 助 于 
认识 气候 变化 对 冰川 和 径流 的 影响 ,而 且 对 冰川 物 
质 平衡 的 预 估 具 有 重要 的 参考 价值 。 目 前 关于 祁 
连 山 冰川 物质 平衡 的 监测 仅 限于 老虎 沟 12 号 冰川 、 
七 一 冰川 八 一 冰川 等 为 数 不 多 的 几 条 冰川 ""。 
此 外 ,冰川 物质 平衡 的 观测 需要 耗费 极 大 的 人 力 物 
73 ,观测 年 限 较 短 且 不 连续 ,观测 时 段 仅 限 于 20 世 
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P 段 和 东 段 流域 的 冰川 消融 速率 大 于 西 段 的 冰川 。 


纪 70 年 代 以 及 2000 年 之 后 ,缺乏 长 时 间 序 列 的 变 
化 研究 。 虽 然 , 目前 在 祁连山 开展 了 一 些 冰 川 物质 
平衡 的 相关 的 工作 ,例如 用 遥感 监测 和 大 地 测量 法 
手段 研究 冰川 面积 和 冰 储 量 的 变化 *。 此 外 ,已 
有 人 研究 也 利用 度 日 模型 重建 祁连山 的 冰川 物质 平 
衡 变化 ,但 是 上 述 研 究 缺 乏 深入 具体 到 祁连山 
流域 尺度 的 冰川 物质 平衡 变化 研究 ,对 气候 一 冰川 
的 关系 认识 仍然 不 足 。 度 日 因子 模型 计算 简单 ,所 
需 参 数 容易 获取 ,在 流域 尺度 上 能 够 得 到 与 能 量 平 
衡 模 型 类 似 的 结果 ”1 。 本 研究 利用 加 入 辐射 项 的 
分 布 式 度 日 模型 ,重建 祁连山 流域 尺度 冰川 的 历史 
物质 平衡 变化 序列 ,从 而 有 利于 我 国 西部 地 区 在 冰 
冻 圈 变化 影响 下 制定 科学 合理 的 水 资源 管理 措施 和 
实施 方案 ,实现 区 域 环境 和 经 济 的 可 持续 发 展 ”。 


1 研究 区 概况 


祁连山 位 于 青藏 高 原 东北 缘 ( 图 1), 东 起 乌 鞘 
岭 , 西 至 当 金 山口 , 南 靠 柴 达 木 倪 地 , 北 临河 西 走 
廊 , 长 约 800 km, 宽 约 300 km” o ZÆ L 4b X JE F 
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图 1 祁连山 区 流域 及 气象 站 点 分 布 示 意图 
Fig.1 Schematic diagram of the distribution of watershed and meteorological stations in the Qilian Mountains 
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影响 , 西 段 受 西风 环流 的 控制 。 冬 半年 主要 受到 蒙 
古 高 压 的 控制 ,天 气 寒冷 且 干 燥 ; 夏 半年 受到 来 自 
印度 洋 的 西南 季风 和 来 自 太平 洋 的 东南 季风 共同 
影响 ,年 均 温 5% ,年 降水 量 为 230 mm ,其 中 降水 主 
要 集中 在 夏季 , 约 占 全 年 降水 量 的 70% ,空间 分 布 
由 东南 向 西北 逐渐 减少 ?0 。 

根据 第 二 次 冰川 编目 统计 的 数据 2 ,祁连山 共 
有 冰川 2683 条 ,面积 (1597.81+70.3)km ,储量 为 
(84.48+3.13 )km ,集中 分 布 在 99°E 以 西海 拔 4800~ 
5200 m 以 上 的 高 山地 带 ( 图 1)。 位 于 下 勒 河流 域 的 
老虎 沟 12 号 冰川 是 祁连山 最 大 的 山谷 冰川 ,面积 ; 
20.4 km ,是 我 国 进行 冰川 监测 研究 的 第 一 个 野外 
观测 站 ,也 是 本 研究 中 优化 参数 的 参考 冰川 。 其 中 
玻 勒 河流 域 的 冰川 数量 最 多 且 面 积 最 大 ,发 育 有 
660 条 , 占 总 数量 的 25% ;面积 为 518 km? , 占 总 冰川 
面积 的 32%。 冰 川 数量 最 少 、 面 积 最 小 的 是 巴 音 郭 
勒 河流 域 ,只 有 10 条 冰川 , 占 总 数量 的 0.4% ;面积 
为 0.18 km , 占 总 冰川 面积 的 0.01%( 图 2)。 


2 研究 方法 和 数据 
2.1 方法 


2.1.1 分 布 式 度 日 模型 ”基于 简单 的 温度 指数 度 日 
模型 (degree-day temperature-index models, DDM ) 通 
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图 2 祁连山 各 流域 冰川 统计 


Fig. 2 Glacier statistics of Qilian Mountains 


过 增加 辐射 项 来 计算 冰川 消融 , 相 比 于 简单 的 温度 
ra BOISE 日 模型 ,增加 辐射 项 的 分 布 式 度 日 模型 能 
较 好 的 体现 冰川 物质 平衡 的 空间 变化 ”: 


MB=P-DDF-PDD+ asof) (1) 


式 中 :MB 是 冰川 年 物质 平衡 (mm w.e.) ; DDF Z& E H 
因子 (mm w.e. «d C20; PDD Æ ERC-a’); P 
是 固态 降水 量 (mm); a I b ESAT IF sr Fa JI] 4: 
AS BAAR HE C BI aT AE AW em); FEV 
表面 所 有 栅 格 单元 格 的 平均 年 短波 辐射 (W.m”)。 

2.1.2 短波 辐射 短波 辐射 模型 方法 参考 Kumar 
等 ,模型 考虑 到 纬度 ,高程 .坡度 . 坡 向 以 及 地 面 
反射 率 等 因素 对 其 的 影响 。 包 括 直接 辐射 .散射 辆 
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射 和 反射 辐射 。 (3) 反射 辐射 
(1) 直接 辐射 反射 辆 射 的 大 小 取决 于 表面 的 倾斜 度 和 地 面 


太阳 辐射 穿 过 大 气 层 时 会 衰减 ,假设 天 空 为 晴 

空 可 以 用 布 格 定律 ”: 
lla" (2) 
式 中 : 五 和 万分 别 是 地 面 和 大 气 上 界 接受 的 辐射 
(W.nm2); 天 是 吸收 常数 ;好 是 空气 质量 比 。 太 阳 高 
度 角 为 90" 时 W=1, 太 阳 在 地 平 线 时 W=30。 
I, = S,[1 + 0.0344 cos(360°N)/365] 


(3) 
M=|1299+ (614 sin ay] -614sina (4) 
式 中 : SV 为 太阳 常数 (1367 Wem’); a 为 太阳 高 
度 角 。 


倾斜 表面 的 太阳 辐射 7: 
1,=1,cosi (5) 
/为 潜在 的 地 表 短波 辐射 ,计算 方法 如 下 : 
1 =1r (6) 


cosi = sin 8 (sin L cos — cos L sinf cos a,) * 
cos6, cosh (cos L cos + sin L sinB cosa, ) + (7) 
cos, sinB sin a, sinh, 
式 中 :i 是 表面 法 线 与 太阳 方向 的 夹 角 ; B 是 坡度 ; 
a, 是 地 表 方 位 角 ; 6, 是 太阳 偏 角 ; 工 是 纬度 ; h 是 
小 时 角 。 


的 反射 系数 ,计算 公式 如 下 ”: 
[,=rl,z, sin B/2sina (12) 
7, =0.217 + 0.706r, (13) 
RP: 了 7 为 反射 辐射 ; r 是 地 面 反射 系数 ; r, 是 反射 
透射 率 。 
2.2 数据 及 处 理 
模型 中 的 冰川 输入 数据 为 国家 冰川 冻 土 沙漠 
科学 数据 中 心 的 中 国 第 二 次 冰川 编目 数据 (http:/ 
www.nedc.ac.cn) "; DEM 数据 为 SRTM 高 程 数据 ， 
分 辨 率 为 900 m。 降 水 数据 来 源 于 “国家 青藏 高 原 科 
学 数据 中 心 ” 的 中 国 高 寒山 区 月 降水 数据 集 
(CAPD) (http://data.tpde.ac.en ) ^**" , 空间 分 辨 率 为 
1000 m。 为 了 方便 计算 ,在 ArcGIS 软件 中 将 其 统一 
重 采 样 为 90 m 的 分 辨 率 并 处 理 为 年 降水 量 。 气 温 
数据 来 源 于 国家 气象 站 和 老虎 沟 12 号 冰川 的 气象 
站 (图 1)。 冰 川 物质 平衡 的 实测 数据 主要 来 自 老 虎 
沟 12 号 冰川 和 七 一 冰川 。 
2.3 模型 的 参数 率 定 
正 积 温 根 据 自 动 气 象 站 的 观测 值 与 海拔 建立 
线性 关系 并 外 推 到 冰川 格 点 上 (图 1)。 考 虑 到 空间 
异 质 性 ,将 祁连山 划分 为 西 段 .中段 和 东 段 ,在 各 段 
分 别 建立 关系 得 到 正 积 温 的 海拔 递减 率 分 别 为 


5, = 23.45 sin]360°(284 + N)/365] (8) i 
c WRN APR AIRE: -1.08 C+m'',-1.02 Tm ! 3-143 Cm", HYHA 
obe 化 上 选取 各 流域 内 或 距离 较 近 的 代表 气象 站 , 即 玉 
uL NEL Rta rs wee tet Spt ce ae 
(2) 散射 辐射 mE 站 1961—2013 年 的 变化 与 各 栅 格 点 建立 关系 ( 表 
散射 辆 出 根 据 Cate 计 算 “”: 1)。 度 日 因子 除 北大 河流 域 根据 七 一 冰川 的 实测 数 
li=hracos B/2sina CO 据 率 定 ,其 他 流域 均 根据 老虎 沟 12 号 冰川 的 实测 数 
7,=0.217 -0.2497, QD 据 建立 与 海拔 的 线性 关系 进行 插值 ( 表 2)。 固 态 降 
式 中 : 7 为 散射 太阳 辐射 ;zi 是 辐射 散射 系数 。 水 占 比 为 0.9”。 辐 射 因 子 在 时 空 上 假设 为 常数 忆 ， 
表 1 国家 气象 站 站 点 信息 
Tab.1 Information of national weather station 
站 点 纬度 (N) AREE) 海拔 /m 流域 
玉 门 40°16'12" 97°1'47" 1527 HAA 
野牛 沟 38°25'12" 99°34'48” 3320 黑河 青海 湖 
Bde 37°12! 102°52'12" 3044 石 羊 河 
KEHE 37°51" 95°22" 12" 3174 APAI EEE ` 鱼 卡 - 塔 塔 楼 、 哈 尔 腾 河 
酒泉 39°46'12" 98°28'48" 1478 北大 河 
门 源 37°22'48" 101°37'12" 2851 大 通 河 
冷 湖 38°45’ 93°19'48" 2770 党 河 
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X2 度 日 因子 . 正 积温 随 海拔 的 递减 率 
Tab. 2 Change of degree-day factor and positive 


accumulated temperature with altitude 


流域 递减 率 
EH ” 疏 勒 河 流域 与 其 他 流域 —0.0293Ele+147.97 
因子 北大 河流 域 -0.0309Ele+147.97 
正 积温 中 有 段 ( 青 海 湖 和 黑河 流域 ) -1.02 Cem" 


东 段 ( 石 羊 河流 域 和 大 通 河流 域 ) -1.15 *C* m" 
西 段 (其 他 流域 ) -1.08 C- m^ 
LE: Ele 为 海拔 (m)。 


根据 在 老虎 沟 12 号 冰川 的 物质 平衡 观测 结果 对 其 
系数 进行 优化 ,a 取 24,5 取 19.8。 
2.4 模型 的 误差 分 析 

模型 的 误差 源 主要 为 输入 数据 和 参数 的 不 确 
定性 。 为 了 计算 模型 结果 的 不 确定 性 ,分 别 对 输入 
数据 即 正 积温 (+20 %C)、 降 水 (+10%) 和 参数 度 日 因 
子 (+1 mm w.e.d!'.°07 1) 的 物质 平衡 敏感 性 进行 分 
析 。 模 拟 结果 的 总 误差 为 三 者 的 误差 之 和 。 


3 结果 分 析 


3.1 冰川 物质 平衡 变化 

模拟 结果 表明 ,以 1995 年 为 界 祁 连 山 冰川 的 物 
质 平衡 变化 可 划分 为 2 个 阶段 ,1961 一 1995 年 的 冰 
川 物 质 平衡 呈正 负 交 替 变化 且 变 化 幅度 较 小 ,平均 
物质 平衡 为 (0.11+0.13 ) m w.e.*a,1996—2013 年 为 
人 负 平 衡 且 变化 幅度 大 ,年 均 物 质 平衡 为 -0.54+0.13 
m we.a (图 3)。 其 中 2010 年 的 冰川 物质 平衡 最 


1996 一 2013 年 
mmu 1961—19954F 


posi 


物质 平衡 /(m w.e.-a) 


图 3 祁连山 及 各 流域 不 同 阶 段 的 年 平均 物质 平衡 


Fig. 3 Annual average massl balance of different stages in 


the Qilian Mountains and basins 
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小 ,为 (-0.72+0.13)m w.e. , 1983 年 的 物质 平衡 最 大 
为 (0.21+0.13 )m w.e.( 图 4)。1961 一 2013 年 祁连山 
冰川 的 累积 物质 平衡 为 (-12.76+4.24)m w.e., 其 中 
1961—1995 年 间 祁 连 山 的 冰川 累积 物质 平衡 为 轻 
微 的 正平 衡 ,1996 一 2013 年 的 累积 物质 平衡 迅速 下 
降 , 冰 川 消融 强烈 (图 5)。 

祁连山 冰川 物质 平衡 在 20 世 纪 60 一 90 年 代 基 
本 维持 稳定 ,主要 原因 是 正 积温 增加 趋势 微弱 且 降 
水 增加 显著 , 即 降 水 量 的 增加 能 够 弥补 正 积 温 上 升 
带 来 的 冰川 消融 ; 正 积 温 上 升 与 降水 量 下 降 共同 作 
用 导致 90 年 代 冰 川 物质 平衡 迅速 下 降 ;2000 年 之 后 
虽然 降水 量 有 所 回升 但 无 法 弥补 正 积温 显著 升 高 
带 来 的 冰川 消融 (图 6)。 

1961 一 2013 年 巴 音 郭 勒 河流 域 . 哈 拉 湖 流域 和 
党 河流 域 的 冰川 物质 平衡 基本 为 负 , 其 他 流域 则 在 
1961—1995 年 为 下 ,1995 一 2013 年 为 负 。 巴 音 郭 勒 
河流 域 1996 一 2013 年 的 冰川 亏损 最 大 ,平均 物质 平 
衡 达 到 了 (-1.23+0.22)m w.e.…a'。 鱼 卡 塔 塔 棱 流 域 
的 冰川 1995 年 之 后 消融 量 最 少 , 平 均 为 (-0.02+ 
0.06)m w.e.-a ([K] 3) ; 

空间 分 布 上 祁连山 东 段 和 中 段 的 冰川 消融 速 
率 大 于 西 段 (图 4a)。 巴 音 郭 勒 河流 域 的 冰川 物质 
平衡 变化 速率 最 快 ,达到 了 -24.9 mm w.e.…a ,是 因 
为 巴 音 郭 勒 河 流域 的 正 积温 升 温 速 率 最 快 ,为 
7.71 Ca ;其 次 为 石 羊 河流 域 和 大 通 河流 域 ,分 别 为 
(-20.9+0.07 )mm w.e..a 和 (-19.3+0.04)mm w.e.-a’'; 
鱼 卡 塔 塔 棱 流 域 的 冰川 消融 速率 最 慢 为 (-3.9+ 
0.06) mm w.e.-a (€ 3 和 图 4a)。 

巴 音 郭 勒 河流 域 . 哈 拉 湖 流域 和 党 河流 域 的 冰 
川 累 积 物质 平衡 恒 为 负 。 其 中 巴 音 郭 勒 河流 域 的 
冰川 累积 物质 平衡 最 低 为 (-35.1+11.6)m w.e., 哈 拉 
湖 流域 和 党 河流 域 的 冰川 累积 物质 平衡 分 别 为 
(-26.8+7.95)m we. 和 (-12.9+4.24)m w.e.。 青 海 湖 
流域 黑河 流域 和 鱼 卡 塔 塔 棱 流 域 的 冰川 累积 物质 
平衡 为 正 , 分 别 为 (5.1+1.38)m w.e.、(1.9+1.48)m w.e. 
和 (3.9+3.28)m w.e.。 原 因 是 青海 湖 流 域 的 冰川 区 
降水 增加 速率 较 快 (2.2 mm * a7) ;黑河 流域 的 冰川 区 
正 积温 上 升 速率 最 慢 (0.64 Ca!) (483) ; f KEE 
棱 流 域 的 冰川 区 分 布 在 高 海拔 区 。 石 羊 河流 域 、 哈 
尔 腾 河流 域 和 北大 河流 域 的 冰川 累积 物质 平衡 在 
20 世 纪 90 年 代 由 正平 衡 转 为 负 平 衡 。 大 通 河流 域 
的 冰川 累积 物质 平衡 于 2009 年 由 正平 衡 转 为 负 平 
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图 4 1961—2013 年 祁连山 冰川 物质 平衡 时 空 变化 


Fig.4 Temporal and spatial changes of glacier mass balance in the Qilian Mountains during 1961-2013 
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图 5 1961 一 2013 年 祁连山 冰川 累积 物质 平衡 变化 
Fig. 5 Change of accumulative mass balance of Qilian 
Mountains glaciers from 1961 to 2013 
衡 (图 5)。 


3.2 模型 的 验证 与 对 比 

将 老虎 沟 12 号 冰川 物质 平衡 模拟 结果 和 花 杆 
实测 数据 及 DEM (Digital Elevation Model) 差 分 结果 
对 比 。 结 果 表 明 ,2011 年 和 2012 年 的 模拟 值 与 实测 
值 相当 ,20 世 纪 70 年 代 的 物质 平衡 较 实 测 数 据 略 有 
低估 (RMSE=0.08 m w.e.)( 图 7a)。DEM 差分 结果 显 
zs 1957—1989 年 老虎 沟 12 号 冰川 物质 平衡 为 
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图 6 1961 一 2013 年 祁连山 冰川 物质 平衡 及 正 积 温和 
降水 的 年 代 际 变化 


Fig.6 Interdecadal variation of mass balance and positive 


cumulative temperature and precipitation in Qilian Mountains 
glaciers from 1961 to 2013 


-0.33 m w.e., 1989—2015 年 间 物 质 平衡 为 -14.36 
m w.e., 与 本 研究 结果 接近 (图 7b )。 

将 七 一 冰川 和 宁 缠 河 3 号 冰川 的 物质 平衡 模拟 
结果 与 实测 值 对 比 发 现 ,模型 总 体 上 能 够 模拟 出 七 
一 冰川 物质 平衡 的 变化 趋势 (RMSE=0.13 m w.e.) , 
在 20 世 纪 70 一 80 年 代 的 模拟 效果 较 好 ,2000 年 和 
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#3 1961 一 2013 年 各 流域 冰川 区 正 积温 \ 降 水 和 物质 
平衡 的 变化 率 
Tab.3 Change rates of positive accumulated 
temperature, precipitation and mass balance in the glacial 


regions of each basin during 1961-2013 


I 降水 变化 率 ” 正 积温 变化 率 ”物质 平衡 变化 率 
流域 


/(mm:a‘') (C+a') /mm w.e.*a') 
mT 2.12 1.76 -17.0 
青海 湖 2.20 1.99 -16.3 
鱼 卡 塔 塔 棱 2.04 4.58 -3.9 
哈 拉 湖 2.09 5.30 -15.5 
黑河 1.89 0.64 -12.1 
巴 音 郭 勒 河 1.88 7.71 -24.9 
党 河 2.19 1.63 -7.0 
北大 河 1.99 3.75 -16.7 
哈 尔 腾 河 2.13 5.68 -122 
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2001 年 的 物质 平衡 模拟 较 实 测 值 略 偏 高 (图 7c)。 
刘 宇 硕 等 外 利用 探 地 雷达 并 结合 影像 资料 计算 了 
宁 缠 河 3 号 冰川 1972 一 2009 年 冰川 物质 平衡 水 当 
量 大 约 为 -7.99 m w.e., 与 本 人 研究 结果 接近 (图 7d)。 
为 了 进一步 验证 模型 的 准确 度 , 将 本 研究 与 
Yao 等 "在 七 一 冰川 1975 一 2010 年 的 累积 冰川 物质 
平衡 重建 结果 进行 比 对 (图 8a) ,可 以 看 出 ,20 世 纪 
90 年 代 之 后 二 者 的 累积 物质 平衡 结果 较为 接近 ,90 
年 代 之 前 本 研究 的 模拟 结果 低 于 已 有 人 研究。 也 就 
是 说 与 Yao 等 "的 研究 结果 相 比 ,90 年 代 之 前 本 研 
究 的 冰川 物质 平衡 模拟 值 偏 小 ,90 年 代 之 后 的 模拟 
{HK Wang °F] FA CRU TS 气温 和 降水 驱动 数 
据 运 用 度 日 模型 重建 了 高 亚洲 的 冰川 物质 平衡 变 
化 ,本 研究 与 其 在 祁连山 的 研究 结果 进行 对 比 , 二 
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图 7 典型 冰川 物质 平衡 模拟 与 实测 对 比 


Fig.7 Comparison of mass balance simulation and measured 
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8 本 研究 与 已 有 研究 的 物质 平衡 重建 结果 对 比 


Fig. 8 Comparison with the results of mass balance reconstruction from existing studies 
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图 9 2000 一 2010 年 本 研究 与 大 地 测量 法 计算 的 物质 平衡 
结果 对 比 
Fig. 9 Comparison of basin-scale mass balance calculated by 
the geodetic method from 2000 to 2010 


者 在 90 年 代 的 变化 过 程 相 一 致 ,但 90 年 代 之 前 本 
研究 的 模拟 结果 相 比 偏 大 ,90 年 代 之 后 的 模拟 结 

偏 小 (图 8b)。 原 因 是 数据 来 源 和 参数 率 定 方 法 的 
不 同 造 成 。 此 外 ,将 本 人 研究 结果 与 大 地 测量 法 计算 
的 2000 一 2010 年 祁连山 冰川 流域 尺度 的 物质 平衡 
变化 进行 对 比 3, 模 型 的 重建 结果 与 各 流域 的 大 地 
测量 法 结果 相 吻 合 , 均 在 其 误差 范围 之 内 (图 9)。 


4 讨论 与 结论 


祁连山 的 冰川 物质 平衡 于 20 世 纪 90 年 代 中 期 
发 生 突变 ,这 与 很 多 研究 结果 一 致 下。 1995 
年 为 界 可 以 划分 为 2 个 阶段 :1961 一 1995 年 冰川 物 


质 平衡 呈正 负 交 替 变化 ,冰川 积累 量 与 消融 量 基本 
持平 ; 1996—2013 年 为 负 平衡 ,平均 物质 平衡 为 
(-0.54+0.13)m w.e- a', 主要 原因 是 正 积 温 的 逐年 
上 升 导致 冰川 消融 加 剧 。 巴 音 郭 勒 河流 域 的 冰川 
物质 平衡 变化 最 快 为 -24.9 mm w.e. a”, ARERI 
流域 的 冰川 消融 速率 最 慢 为 -3.9 mm w.e.*a '。 

1961—2013 年 祁连山 冰川 的 累积 物质 平衡 为 
-12.76 m w.e.。 其 中 巴 音 郭 勒 河流 域 的 冰川 累积 
质 平 衡 最 低 为 -35.1 m w.e.; 青 海 湖 流域 .黑河 流域 
和 鱼 卡 塔 塔 棱 流 域 的 冰川 累积 物质 平衡 为 正 , 小 型 
冰川 分 布 较 多 的 巴 音 郭 勒 河流 域 . 哈 拉 湖 流域 和 党 
河流 域 的 累积 物质 平衡 则 一 直 为 负 。 

度 日 因子 作为 模型 中 的 重要 参数 ,前 人 的 研究 
多 数 将 度 日 因子 划分 为 十 和 雪 并 分 别 设 定 为 常 
数 。 本 文 的 度 日 因子 是 根据 典型 冰川 的 实测 数据 
与 海拔 建立 线性 关系 ,但 受 限 于 冰川 实测 数据 的 不 
AE ,存在 着 代表 性 不 足 的 问题 ,未 来 随 着 实测 资料 
的 增加 需要 进一步 进行 验证 和 细 化 。 此 外 ,模型 的 
模拟 过 程 中 没有 考虑 到 度 日 因子 随时 间 的 变化 , 例 
如 近年 来 冰川 杂质 的 增多 5 会 导致 度 日 因子 增 大 ， 
因此 模型 会 低估 冰川 的 物质 平衡 量 。 

研究 表明 祁连山 第 一 次 冰川 编目 (1956 一 1983 
年 ) 到 第 二 次 冰川 编目 (2005 一 2010 年 ) 近 50 a 间 冰 
川 面积 变化 为 -20.88%'1。 本 研究 中 物质 平衡 的 冰 
川 输入 数据 为 第 二 次 冰川 编目 的 数据 ,没有 考虑 到 
冰川 边界 及 高 程 变 化 ,导致 模 型 高 佑 冰川 的 物质 平 
衡量 。 但 是 经 过 与 实测 数据 及 其 他 的 研究 结果 对 
比 认 为 模型 总 体 上 能 够 重 现 各 流域 冰川 的 物质 平衡 
变化 ,可 以 为 冰川 水 文 模型 的 模拟 提供 前 提 条 件 。 
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Reconstruction of the glacier mass balance in the Qilian Mountains 
from 1961 to 2013 


WANG Lihui", QIN Xiang’, CHEN Jizu, ZHANG Dongwei’, 
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(1. State Key Laboratory of Cryospheric Science, Northwest Institute of Eco-Enivironment and Resources, CAS, 


Gansu, Lanzhou 730000, Gansu, China; 2. University of Chinese Academy of Sciences, Beijing 100049) 


Abstract: Glacier mass balance is the most direct and reliable indicator of climate change. Its dynamic change 
likely alters the ice storage and runoff of glaciers. Most glaciers in the world have retreated because of global 
warming, and this phenomenon has been accelerating for the last 20 years. The glacier area in the Qilian 
Mountains is the ^wet island" of the arid and semiarid regions of the Qaidam Basin and the Hexi Corridor. In 
these parts, the main replenishing areas of water resources are abundant glacier and precipitation resources, which 
provide abundant and stable water resources for the economic development of these regions. Glacier melt water 
accounts for 25%-60% of the river runoff in the inland river basins of northwest China. As such, changes in mass 
balance in the Qilian Mountains should be monitored and estimated to develop scientific and reasonable measures 
and implementation programs under the influence of cryospheric changes in western China and to establish risk 
control and policy regulations. In this study, a distributed degree-day model with radiation terms was used to 
reconstruct the glacial mass balance from 1961 to 2013 at the watershed scale of the Qilian Mountains. Input data 
included the second glacial catalog data set of China, precipitation data from the National Qinghai- Tibet Plateau 
Science Data Center, and temperature data from weather stations and digital elevation models. The degree-day 
factor was obtained from existing monitoring report on glaciers. Changes in the glacier mass balance in the Qilian 
Mountains could be divided into two stages in 1995. The glacier mass balance exhibited alternate changes 
between positive and negative values, and such variations were relatively small from 1961 to 1995. The amounts 
of glacier melting and accumulation were basically the same. Glacier melting intensified rapidly from 1996 to 
2013. The mass balance remained negative, and the average mass balance was (—0.54+0.13)m w.e.-a ! mainly 
because the annual increase in the positive accumulated temperature led to the accelerated melting of glaciers. 
The cumulative mass balance of glaciers in the Qilian Mountains from 1961 to 2013 was (—12.76+4.24)m w.e., 
the glacier melting rates in the middle and eastern parts of the Qilian Mountains were higher than that in the 
western part. 


Keywords: Qilian Mountains; watershed; glacier mass balance; reconstruction 


